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We study how the properties of the four neutrahno states, Xi (i — li 2, 3, 4), can be investigated 
04 , at the Large Hadron Colhder (LHC), in the case when the lightest one, xi, has a mass < 50 GeV 

and is stable. This situation arises naturally in supersymmetric models where gaugino masses are 
not unified at a Grand Unified (GUT) scale and R-parity is conserved. The main features of these 
f~i ■ neutralino states are established by analytical and numerical analyses, and two scenarios are singled 

out on the basis of the cosmological properties required for the relic neutralinos. Signals expected 
at LHC are discussed through the main chain processes started by a squark, produced in the initial 
' proton-proton scattering. We motivate the selection of some convenient benchmarks, in the light 

,.1^ , of the spectroscopical properties (mass spectrum and transitions) of the four neutralino states. 

Branching ratios and the expected total number of events are derived in the various benchmarks, 
and their relevance for experimental determination of neutralino properties is finally discussed. 
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■ I. INTRODUCTION 
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Light neutralinos, i.e. neutralinos with a mass < 50 GeV, arise naturally in supersymmetric models where 



o 

00 

] gaugino masses are not unified at a Grand Unifi.ed (GUT) scale. When R-parity conservation is assumed, these 
neutralinos offer a very rich phenomenology under various cosmological and astrophysical aspects [l|, Q . 

In the present paper wc study how these light neutralinos can be investigated at the Large Hadron Collider (LHC), 
' which will soon start being operated at CERN [^, Ql . 

To this purpose we first delineate what are the main properties of the four neutralino states, Xi (i = 1; 2, 3, 4), 
in the case when the lightest one, xi (or X iii short), has a mass < 50 GeV. This neutralino is stable, in case 
it occurs to be the Lightest Supersymmetric Particle and R-parity is conserved. In particular, special asymptotic 
schemes for the four neutralino states are analysed: two hierarchical schemes (with so-called normal and inverted 
hierarchy, respectively) and an almost degenerate one. 

Once the properties of all the neutralino states are established by analytical and numerical analyses, two main 
scenarios are singled out on the basis of the cosmological properties required for the relic neutralinos. 

Then, we analyse the signals expected at LHC, in terms of the main chain processes which are initiated by a squark 
and lead finally to xi through an intermediate production of Xi (i = 2, 3, 4). We evaluate the branching ratios for the 
various processes and discuss their features in terms of the spectroscopic properties (mass spectrum and transitions) 
of the four neutralino states. Various benchmarks of special physical interest arc considered. 

The total number of events expected for these processes at LHC is estimated for a typical representative value of the 
integrated luminosity and their relevance for experimental determination of neutralino properties is finally discussed. 

The scheme of our paper is the following. The features of the employed supersymmetric model are presented in 
Sect, ini where also the main structural properties of the four neutralino states are described. In Sect. Illll we delineate 
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the scenarios imposed by cosmology on light relic neutralinos. In Sect. IIVI we derive branching ratios and the total 
number of events expected at LHC. Finally, conclusions are drawn in Sect. |Vl 



II. THE SUPERSYMMETRIC MODEL 



The supersymmetric scheme we employ in the present paper is the one described in Ref. [ij: an 
effective MSSM scheme (effMSSM) at the electroweak scale, with the following independent parameters: 
Ml, M2, Ms, fi,tanP,mA,mq,irn and A. Notations are as follows: Mi, M2 and M3 are the U(l), SU(2) and SU(3) 
gaugino masses (these parameters are taken here to be positive), fi is the Higgs mixing mass parameter, tan/? the 
ratio of the two Higgs v.e.v.'s, niA the mass of the CP-odd neutral Higgs boson, rriq is a squark soft-mass com- 
mon to all squarks, mj is a slepton soft-mass common to all sleptons, and A is a common dimensionless trilinear 
parameter for the third family, ~ Af = Arriq and Af = Am^ (the trilinear parameters for the other families 
being set equal to zero). In our model, no gaugino mass unification at a Grand Unified (GUT) scale is assumed. 
The following experimental constraints are imposed: accelerators data on supersymmetric and Higgs boson searches 
(CERN e+e~ collider LEP2 and Colfider Detectors DO and CDF at Fermilab [^); measurements of the 6^5 + 7 
decay process 0: 2.89 < B{b — > s + ■ 10^ < 4.21 is employed here (this interval is larger by 25% with respect 
to the experimental determination Q in order to take into account theoretical uncertainties in the supersymmetric 
(SUSY) contributions Q to the branching ratio of the process (for the Standard Model calculation, we employ the 

M^) we take 



recent NNLO results from Ref. [y|)); the upper bound on the branching ratio BR{B^ —t fi 
BR{B'^ ji" 4- /i"*") < 1.2 • 10~^; measurements of the muon anomalous magnetic moment 



(.g^ - 2)/2: for the 



deviation Aa^ of the experimental world average from the theoretical evaluation within the Standard Model we use 



here the range —98 < Aa^ • 10^^ < 565, derived from the latest experimental [ll| and theoretical [121 data 



A. Composition of the neutralino states 



The linear superpositions of bino B, wino W'^'^^ and of the two Higgsino states H^, H2 which define the four 
neutralino states are written in the following way: 



X^ = afB + a(*%(3) ^ ^W^o ^ ^yv jjo ^ 2, 3, 4). 



(1) 



These states diagonalize the mass matrix 



Ml 






M2 



-mzsecp mzcecp 
V mzsgsp -mzcgsp 



-mzsgcis mzsesii \ 
mzC0Cf3 -mzcesp 
-Ai 

/ 



(2) 



where sp = sin/3, cp = cos (3, and sg = sinO, cg = cos 9, 6 being the Weinberg angle. The mass eigenvalues (with 
signs) will be denoted by m^. The smallest mass eigenvalue \mi\ will also be denoted by . 
From the set of equations 



{Ml ~ ■mi)a-^'' ~ mzsecpa}^' + mzsgSfjal'' = 
{M2 - mi)a2^ + mzcecpa'^^ - mzcespaf^ = 
—razsgcpa\ + mzCgCfja^ ~ rriia^ — fia\ = u 
mzsgspal — mzcgspa2 — iJ-a^ — mia\ = 0, 

that follow from the diagonalization of the mass matrix of Eq. ([2]) , one obtains the ratios 



(3) 
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Ml - rrii 
M2 - TO,; 



cote 
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ay 



flW ^cp + miSfj 

which, together with the normaUzation conditions, provide the compositions of the four neutrahno states [l^l- 

The LEP lower Hmit on the chargino mass (to^± > 100 GeV) sets a lower bound on both and M2: |/i|,M2 > 
100 GeV. Since, on the contrary, Mi is unbound, the lowest value for occurs when 

~ A/i « 1^1, Afa. (5) 

Thus, X = Xi is mainly a Bino, whose mixings with the other interaction eigenstates are readily derived from Eqs. 
(HD to be 

(1) c (1) (1) 

-—cote, -rjTC^SgSf} , —J--— ■ , (6) 

where ^1 = mi — Mi. Notice that, in deriving Eq. ([5]), we have also taken into account that in the present paper we 
will only consider tan/3 > 10. 

From Eqs. (|4]) one also obtains approximate expressions for the compositions of the eigenstates which correspond 
to the asymptotic mass eigenvalues: — ±fi and mi — A/2. That is: 

a) for the neutralino states Xi with m,; ~ ±/i, 

^ ±/x a^^ 26ge(±/i-M2) 4'> £2 
cote, -7-7 - TTTT-TT-^Tr^TTT' -777-^1 + —, (7) 



where £2 = ±M ~ '^i- 

b) for the neutralino state Xi with mi ~ A/2, 

af A/2 '''''' 4) A/z(M2C/3+/iS/3)c2A/2' 4) A/2S0 + MC^' ^' 

where £3 = A/2 — mi . 

For illustrative purposes, in Figs[lH2]wc show the results of the numerical diagonalization of the matrix of Eq. ([2]) 
versus A/2, for the representative point defined by the following values of the other supcrsymmctric parameters which 
enter in the mass matrix of Eq. ([2]): A/i = 25 GeV, n = 300 GeV, tan/3 = 10. Fig. [1] displays the behaviour of the 
mass eigenvalues |mi|. Fig. [2] the compositions of the states Xi- 

The approximate analytic expressions in Eqs. (I6][8]) and the numerical results in Fig. [2] display quantitatively the 
following properties: (i) xi is mainly a B-ino whose mixing with is sizable at small /z, (ii) X3 h^^s a mass \m^\ ~ 
with a large — H2 mixing, independently of A/2, (iii) X2 and X4 interchange their main structures depending on 
the value of the ratio |/i|/A/2: X2 is dominantly a W-ino (with a sizable subdominance of Hi) for A/2 << |/.t| and a 
maximal HI — H2 admixture for A/2 >> \fi\, whereas Xi is a maximal H° — H2 admixture for A/2 << and a very 
pure W-ino for M2 >> \fJ.\- Also their properties relevant to the case A/2 ^ are transparent from Eqs. (OH]) and 
Fig.l 

Depending on the relative values of the parameters A/2 and fi, it is useful, for the discussion to be developed later, 
to define the following neutralino spectroscopic schemes (notice that we always assume (A/i << A/2, \fi\): (i) normal 
hierarchical scheme (A/2 < \fi\), (ii) degenerate scheme (A/2 ~ (iii) inverted hierarchical scheme (A/2 > \fi\). 

These schemes are depicted in Fig. [3] 
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FIG. 1: Neutralino masses as functions of M2 for Mi = 25 GeV, /i = 300 GeV and tan/3 =10. 
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FIG. 2: Compositions of the neutralino states Xi a-s functions of A/2 for Mi = 
left panel: xii up-right panel: X2, bottom-left panel: xsi bottom-right panel: 
short-dashed lines I03I, long-dashed lines I04I. . 



25 GeV, n = 300 GeV and tan/3 =10. Up 
X4. Solid lines denote |ai|, dotted lines \a2 
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FIG. 3; Asymptotic spectroscopic schemes for the x-i (* = 15 2,3,4) neutralino states . 



B. Cosmological properties 



As discussed in Ref. [l[ , under the assumption that R-parity is conserved, the lower bound on the neutrahno mass 

is provided by the upper bound on the rchc abundance for cold dark matter (CDM), Q,cDMh? ■ 
We recall that the neutralino relic abundance is given by 



3.3 • 10-38 cm2 



^xh'= 7 , (9) 



where < (JannV > = ^''/(o'ann v)int, (cann ^^)int being the integral from the present temperature up to the freeze-out 
temperature Tj of the thermally averaged product of the annihilation cross-section times the relative velocity of a pair 
of ncutralinos, Xf is defined as cc^^ = ^ and g-i,(xf) denotes the relativistic degrees of freedom of the thermodynamic 

bath at Xf. For (iJann v) we use here the standard expansion in S and P waves: ((Tann w) — a -I- b{2xf)^^. 

A host of cosmological observations imply that the cold dark matter content has to stay in the range 0.092 < 
^CDAih^ < 0.124 [1J|. Thus the supersymmetric configurations have to satisfy the cosmological constraint Q^h^ < 

{^CDMh'^)max = 0.124. 

III. SCENARIOS FOR LIGHT NEUTRALINOS 

In Ref. [l| it is shown that the cosmological condition fl^h^ < {flcDMh^)max ~ 0.124 provides a lower limit > 

7 GeV. Such a situation occurs when Mi ~ 10 GeV and ((Jann w) — a receives a sizable contribution by the exchange 
of the A Higgs boson in the s channel. This, in tmm, happens when the two following conditions are met: a) tua is as 
small as its experimental lower bound, m.A —90 GeV, b) the B-ino component of the xi configuration is maximally 

mixed with the H? component (i. e. -jyy ~ 0.4). From the second expression in Eq. ^ one sees that condition (b) 

is satisfied when is small (|^| ~ 100-200 GeV). Moreover, it turns out that tan/3 must be large (tan/3 30 - 45) in 
order for the annihilation cross section through the exchange of the A Higgs boson to be sizable. Finally, the trilinear 
coupling is mildly constrained to stay in the interval — 1 < A < +1. 

In the following we will characterize such a set of supersymmetric parameters as Scenario A. More specifically, 
this scenario is identified by the following sector of the supersymmetric parameter space: Mi ^ 10 GeV, |^| ^ (100 - 
200) GeV, niA ^ 90 GeV, tan/3 '--^30 - 45,-1<j4< +1; the other supersymmetric parameters are not a priori fixed. 

From Eqs. ^ it turns out that, in this scenario, the following hierarchy holds for the coefficients a|^'' of xi 

|af)|>|a«|»|aW|,|a«|. (10) 

It is proved in Ref. [H that at small mx, when > (200 - 300) GeV, the cosmological constraint ft^h^ < 
{^CDMh'^)max IS Satisfied because of stau-exchange contributions (in the t, u channels) to < aannV >, provided that: 
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scenario 


Ml [GeV] 


ImI [GeV] 


tan P 


mA [GeV] 


m-, [GeV] 


A 


~ 10 


100-200 


30-45 


~ 90 




B 


~ 25 


> 500 


< 20 


> 200 


100-200 



TABLE I: Scenarios for light neutralinos as described in Section Hill In scenario A the slepton soft mass mi is unconstrained: 
in our analysis a few representative values for mj- are considered within its natural range 115 GeV < m,- < 1 TeV. Moreover, 
in scenario A: -1 < A < +1, in scenario B: -2 < ^ < +2. 



(i) TOf is sufficiently light, mf ^ 90 GcV (notice that the current experimental limit is mf ^ 87 GeV) and (ii) xi is 
a very pure B-ino {i.e. (1 — a^"'^') ~ 0(10^'^). In such a situation, the lower limit on the neutralino mass is mx > (15 
- 18) Gev 0. 

Let us first discuss the implications of the prerequisite (i). The experimental lower bounds on the sneutrino mass 
and on the charged slepton masses of the first two families imply a lower bound on the soft slepton mass: > 115 
GeV. In order to make the request rrif ~ 90 GeV compatible with mj > 115 GeV, it is necessary that the off-diagonal 
terms of the sleptonic mass matrix in the eigenstate basis, which are proportional to /itan/3, are large. Numerically, 
one finds tan/3 ~ 5000 GeV. 

On the other side, the condition (ii) requires that Cg^'/'^'^'' ^5 10^^, «• e., according to the second expression of 

^(1) 

Eq. ([6]), -jfj ~ SQSp^^ < 10^^. Combining this last expression with the condition |/i|tau/3 ^ 5000 GeV, one finds 

that and tan/3 are bounded by: > 500 GeV, tan/3 < 10. These bounds are somewhat weaker for values of the 
neutralino mass larger than ~ 15-18 GeV. 

The previous arguments lead us to introduce a new scenario, denoted as Scenario B, identified by the following 
sector of the supersymmetric parameter space: Mi ^ 25 GeV, \fi\ > 500 GeV, tan/3 < 20; mj > (100 - 200) GeV, 
—2.5 < A < +2.5; the other supersymmetric parameters are not a priori fixed. 

From Eqs. ^ it turns out that, in this scenario, the following hierarchy holds for the coefficients al^-* of xi 



(11) 



The features of scenarios A and B are summarized in Table U] 



IV. SIGNALS AT LHC 

If kinematically accessible, squarks and gluinos are expected to be copiously produced in the pp scattering process 
at the LHC. In the present paper we limit our discussion to the case in which the gluino is heavier than the squark; 
for definiteness, we set the SU(3) gaugino mass at the representative value M3 = 2 TeV and the squark soft-mass at 
the value niq = 1 TeV. Notice that these two parameters are irrelevant in the specification of the scenarios previously 
defined. 

Many experimental studies are available showing that the production of 1 TeV squarks can be easily discovered at 
the LHC through inclusive analyses based on the request of a high multiplicity of hard jets ant E^^'^", see e.g Q and 



In order to check the supersymmetric model described in the previous sections, the SUSY parameters need to be 
measured. For this purpose, we employ here a strategy which has been developed for the measurement of the masses 
of the SUSY particles based on the identification of exclusive decay chains consisting in sequences of two-body decays 
0j 0, 12 , 13 ■ The most promising decay chains considered in the literature are: 



iff qffXi (sequential chain), 



(12) 



and 



q{Z, h, H, A)xi ^ qf fXi (branched chain). 



(13) 
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where / stands for a fermion, and the neutrahno subscript i can take the values 2, 3 or 4. From the experimental point 
of view, given the large multiplicity of QCD jets present in pp events, the only interesting decays which can be used for 
the identification of exclusive chains are the ones involving light charged leptons (e and /x) , the hadronic decays of r 
(r-jets) and the fragmentation of b quarks (6-jets), which can be experimentally separated from the background. The 
sequential chain with two e or in the final state is particularly useful, as the whole chain consists of three successive 
2-body decays which can be measured very well, and provides enough constraints to allow a model-independent mass 
determination. The sequential chain with two r-jets allows a much less clean measurement, as neutrinos are produced 
in the decay of the r which make the event kinematic less clear, and a much higher jet background is present. It 
is however important to detect this chain as well, because by comparing its rate with the rate for the e, ^ chain, 



information can be extracted on the f mixing 19| . Finally, the branched chain provides less constraints on the SUSY 



masses, as the invariant mass of the two final-state fermions just shows a peak at the value of the resonance appearing 
in the decay. In this case the experimentally interesting decays are the ones in b pairs, especially for the Higgs bosons, 
and in e and /i for the Z. 

These two decay chains have in common the first step consisting in the squark decay q qxi ■ Two types of couplings 
occur in this process: (a) gauge couplings, universal in the quark flavor, and (b) Yukawa couplings, hierarchical in the 
quark mass. Gauge (Yukawa) couplings are proportional to the gaugino (Higgsino) components in x^, respectively. 
The actual number of the qxi states produced in a pp collision depends on the product of the cross section for the 
production process u^pp qq^ qq* , gg, qg) times the branching ratio BR(g qXi)- In this product the contribution 
of the heavy quarks is strongly suppressed, because of their scarcity in the proton composition. Moreover, since the 
relative importance of the Yukawa couplings as compared to the gauge couplings depends on the ratio mg/mz {rriq 
and mz being the quark mass and the Z-boson mass, respectively) the production of Xi is more sizable when gauge 
couplings are effective, i. e. when Xi has large gaugino components. To account for these properties it is convenient to 
define an effective branching ratio BR(g — > qXi), i- e- an average of the branching rations BR(^ 9Xi)'s over the light 



quarks 20|, including both qn and q^. To be conservative, we will use as an effective branching ratio BR((7 qXi) 



the average over the four lightest quarks. From now on we will simply denote this average as BR(g qXi)- 

We turn now to a discussion of the decay process for Xii which takes different routes in sequential and branched 
chains, respectively. In the sequential case, since we have taken niq = 1 TcV, the decay can only proceed through 
a slepton: Xi ^ ^ ^ l^Xi with a branching ratio BK{xi II Hxi) = BR(xi H) BR(Z Ixi)- The size of 
BR(xi — > U) depends sensitively on the Xi composition. If Xi is dominantly a gaugino, because of the universality 
of the gaugino couplings, the branching ratios BR(xi — > H) for the three lepton flavours are about the same ; if Xi 
is dominantly a Higgsino, Xi decays predominantly into a f r pair. In the branched chain Xi decays either through 
the Z-boson or through a Higgs boson. The first case, i. e. Xi ~^ Z + xi, involves only the Higgsino components 
of the two neutralino states; the Z boson subsequently decays into all (kinematically possible) // pairs according to 
the Standard Model branching fractions. The second case, i. e. Xi ~^ {h, A, H) + xii in order to have a sizable BR, 
requires that one neutralino state is dominantly a gaugino, the other dominantly a Higgsino. Since in our scenarios 
Xi is dominantly a B-ino state, Xi ~^ (^i H) + xi is of interest when Xi is dominated by the Higgsino components. 
Because of the hierarchical character of the Yukawa coupling, the subsequent decays of the Higgs bosons arc dominated 
by the production of a 6 - 6 pair. We note that in the case of a pronounced hierarchical (inverted) scheme (see Fig. 
[3]) the direct X4 ~^ (h, H, A, Z)xi decay is suppressed because both Xi and xi are dominantly gauginos. This entails 
that the "long" chains X4 — * {h, H, A, Z)x2,3 ^ {h, H, A, Z)xi can give a sizable contribution. Indeed, for both the 
sequential and branched decays, longer multi-step decay chains, of the type e.g. Xi ~^ ^X2,3 XX' xi, are in 
principle very interesting as they provide additional kinematic constraints. In practice, the superposition of many 
decays with the same final state may be extremely difficult to disentangle experimentally, and will thus add confusion 
rather than information. For this study we will therefore limit ourselves to studying the branching ration for direct 
decays. 

In the following sections, for the sequential decays we will consider the two cases / = e, r. We will first analyze e, 
for which the two chirality states and will be considered separately. The decay Xi ~* ce is largely dominated 
by the gaugino components in Xi] thus, if the Xi decays into both e"^ and are kinematically allowed, the ratio 
of the two branching ratios scales approximately as BR{xi ^ <^^l) I BR{xi ee^) ~ 1/4(1 -f cotga^^ /a[''')'^ , where 
/ can be approximated by Eqs. (|7l8p . This implies that Xi decays dominantly to when it is a pure wino 
(flj^V^i*^ ^ 1); while either of the two selectron final states e^^^ may be important for other compositions. For the 
T case we will consider both the mixed (cigcnmass) states fi and T2. The presence of both states in the decay can 
be hard to disentangle, and thus make the mass and rate measurements more difficult. For the branched decays we 



FIG. 4: Qualitative schemes for benchmarks in sequential decay chains: ^-seql, ^-seq2, B-seq. For each benchmark, extremes 
values for M2 are considered: M2 ~ and M2 > \fi\ for scenario A; M2 < \fi\, M2 ~ and M2 > for scenario B. . 



benchmark 


Ml [GcV] 


/i [GeV] 


tan/3 


ruA [GeV] 


[GeV] 


^-seql 


10 




110 


35 


90 


150,300,500,700 


^-seq2 


10 




150 


35 


90 


120 


^"brc 


10 




110 


35 


90 


150 


i3-seq 


25 




-500 


10 


1000 


120 


Z3-brcl 


25 




-500 


10 


200 


120 


Z3-brc2 


25 




-500 


10 


1000 


120 



TABLE II: Benchmarks for light neutralinos, singled out within the two scenarios A and B of Table U All resuhs will be 
presented as functions of M2 , M2 being varied in the range 100 GeV< M2 < 1000 GeV. The other supersymmetric parameters 
are set at the representative values: Ms — 2 TeV, rriq = 1 TeV, A = 0. 



consider the decay into b pairs for the Higgs bosons and the decays into b and e pairs for the Z. In most of the 
considered models the Z and the SUSY Higgscs are almost degenerate, and, given the experimental resolution on the 
bb peak, cannot be separated. It is therefore very useful to have also the decay into leptons, which will allow the 
experiments to determine the presence of SUSY Higgses in the decay chains. 

Based on this discussion, we report our numerical results for all the decay branching ratios relevant for the decay 
chains of Eqs. [T5|) in a number of benchmarks within the scenarios A and B, previously defined (sec Table ITT)) . 
In these benchmarks, all values of the supersymmetric parameters, except M2, are fixed. The branching ratios, given 
as functions of M2, have been calculated by using the ISASUSY code 2l|, 



A. Sequential chain benchmarks 

In scenario A the parameter \fi\ is required to be /i ^ (100-200) GeV, while slepton masses are unconstrained. On 
the other hand the sequential chain is sensitive to the hierarchy between and toj, since, when > the decay 
of the X2,3 states is not allowed by kinematics, and the process (|12p can proceed only through X4- On the contrary, 
when mj < the three states Xi (* = 2,3,4) take part in the chain 

We then fix the following benchmarks for scenario A: 



^-scql : ^ = 110 GeV = 150, 300, 500, 700 GeV (14) 
yl-scq2 : fi = 150 GeV = 120 GeV, (15) 

where, in both cases, Mi=10 GeV, tan /3=35, niA =90 GeV and A=Q. These two benchmarks are depicted qualitatively 
in FigUand summarized in Table HIl 

For the benchmark ^-seql the branching ratios for the sequential decay chain arc shown in Figs l5l6l where, as 
already mentioned, only the decay of Xi is kinematically allowed. Each panel of the figure corresponds to a different 
value of nil among those of Ea. p4)) : mj =150, 300,500,700 GeV. The notations for the various curves are explained 
in the figure caption. 

The main features of Figsl5]and[6]are readily understood. The Branching Ratio (BR) for q qxi is ~ 15%, since 
the X4 is dominantly wino, and therefore the qn has no BR into it, and the qL decays 60% into a chargino X2 and 
30% into X4, according to the left-handed couplings of the wino. In the regime M2 » mj, the Xi decays with 40% 
BR into W, Z and Higgses because of its non-zero Higgsino component, 60% BR into the left-handed sleptons. So 
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BR(q^q e e xj (sequential, benchmark ^ — seql) 

500 1000 500 1000 

1 



mT=150 GeV 



0.01 



03 

a 
o 

a 

g 0.1 
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o.oo: 
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FIG. 5: Branching ratios for the sequential process q eexi in benchmark ^-seql as functions of Af2- Each panel corresponds 
to a different value of mj. The dashed lines show the branching ratio for the process q —> qxi\ the thin-dotted lines denote 
the branching ratio for the process X4 ~^ e^i ~^ sexi', the thick-dotted line in the top-left panel denotes the branching ratio 
for the process X4 — > ecu — > eexi (the corresponding curves cannot be seen in the other panels because they are too low); the 
thick solid lines denote the branching ratio for the whole sequential decay chain q — > qxi — > gee — > geexi • 



for cl the BR is ^10%. For each of the two f states, which are in this case rather similar in composition and mass, 
the BR is ^ 5% for each state. For instance, for mj ~ 150 GeV: to^j = 133 GeV, uif^ = 176 GeV and for mj ~ 300 
GeV: m^i = 192 GeV, m?, = 314 GeV. 

The different quadrants of Figs. [5] and [S] show by how much the branching ratios arc suppressed for growing values 
of mj, due to the reduction of the available phase space for the decay X4 ^ On the other hand, a growing value of 
nil does not affect BR{1 Ixi)^ which is very close to 1 in all cases. Note that in this benchmark whenever M2 > |^| 
one has m^^ ~ M2, thus a measurement of m^^ would typically give direct access to the value of M2. 

The phenomenology of the sequential decay becomes richer in the benchmark ^-seq2, which is shown in Figs[7]and 
[51 In fact in this case the process (fT2)) can proceed through the production and decay of any of the Xi states. The 
corresponding separate contributions to the decay branching ratios for i = 2,3,4 are shown in each of the first three 
panels of Figsl7]and[8]while the last panel shows the combined total branching ratios. 

The branching ratios for the process q — > qXi, shown in the bottom-left panel of Figs[7H51 are obviously very similar 
the ones already displayed in FigslSHSl and previously commented. As compared to these, the branching ratios for the 
processes q qXi (« = 2, 3) are somewhat suppressed, due to the prevalent Higgsino compositions of X2, X3- For the 
same reason, the BR of X2 and X3 into e^^ji, is at most a couple of percent. The contrary happens for the branching 
ratios BR(xi — > ttxi). In this case X2,X3 decay into fr with larger BRs than X4- because the dominance of the 
Higgsino components in X2j Xs favors their decays into third generation leptons. Notice that in the case of X27 only the 
lighter stau contributes, since f2 is heavier than X2; in the case of xs, the f2 contribution is very strongly suppressed as 
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FIG. 6: Branching ratios for the sequential process q frxi in benchmark ^-seql as functions of Ad2- Each panel corresponds 
to a different value of mj. The dashed lines show the branching ratio for the process q —> qX4', the thin-dotted lines denote the 
branching ratio for the process X4 —> tti ttXi i the thick-dotted lines denote the branching ratio for the process X4 ^ ''"'^2 —> 
TTXi', the thin-solid lines denote the branching ratio for the whole sequential decay chain q qxi — > qrfi qrfxi', the solid 
lines with intermediate thickness denote the branching ratio for the whole sequential decay chain q qx4 qrn qrfxi', 
the thickest solid lines denote the total branching ratio for the whole sequential decay chain q qx4 ^ qrf qTfxi ■ 



compared to the fi contribution again for phase-space reasons. All in all, in the full process q qxi qf Qtyxi, 
the most relevant contribution turns out to be provided by X2- 

The scenario B is characterized by a heavy pseudoscalar Higgs mass (to^ > 200 GeV) and by a very light stau, which 
is required in order to keep the ncutralino relic abundance below the observational limit. Moreover, as explained in 
Section Hi B1 must be large. We define in this case the following benchmark, corresponding to the lightest possible 
mass of the neutralino xi- 

fi-seq : Afi = 25 GeV ^ = -500 GeV tan/? = 10 

TOf = 87 GeV {mi = 120 GeV) A = 0. (16) 

This benchmark is depicted qualitatively in Fig|4]and summarized in Table HIl The branching ratios for this benchmark 
are shown in Figs. [5HTU1 In this case the possible hierarchy between M2 and |/^| is richer than in the previous benchmark 
yl-seql, since now also M2 << can occur. In particular, this implies that the compositions of X2 and Xi Aip the 
one into the other in going from M2 < to M2 > (an example of this behaviour can be seen in FiglH although 
for a slightly diff'erent set of supersymmetric parameters). Here, for M2 < 500 GeV, X2 is mainly a gaugino, whereas 
X4 is dominantly a Higgsino; the other way around, for M2 > 500 GeV. While X2 and xa exchange their roles as M2 
runs over its range, ^3 is steadily a Higgsino. independently of il/2. On the basis of these properties one understands 
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BR(q^q e e xj (sequential, benchmark ^ — seq2) 
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FIG. 7: Branching ratios for the sequential process q ee^i in benchmark ^-seq2 as functions of M2. Each of the first 
three panels corresponds to a different intermediate neutralino state (Xii^ = 2,3,4) in the sequential chain, while the last 
panel refers to combined branching ratios. The dashed lines show the branching ratio for the process q —> qXA\ the thin- 
dotted lines denote the branching ratio for the process X4 eeh — > eexi; the thick-dotted lines denote the branching ratio 
for the process X4 een — > eexi The thin-solid lines denote the branching ratio for the whole sequential decay chain 
q qxi qeeL — > qeexi', the solid-lines with intermediate thickness denote the branching ratio for the whole sequential decay 
chain q qx4 — > qeeji qeexi', the thickest solid lines denote the total branching ratio for the whole sequential decay chain 
9 IXi ~* ~^ ^eexi (some curves do not appear in all panels because they are too low) . 

the features of the branching ratios for the production of the Xi intermediate states (dashed lines). It is also clear 
why their combined results (dashed line in the bottom-right panel) have a much milder dependence on M2. 

As for the peculiar behaviour of the branching ratios for the decays Xi ~^ ^^Xi (dotted lines), notice that their 
sudden drop at low il/2 in the cases of X3 and Xi is due to the opening of some competing decay channels. In fact, 
for AI2 < \^J-\ one has ~ 7714 ~ \fi\ and, at the same time, the chargino mass is of order M2. Thus, under these 
circumstances, xs s-^d X4 have a sizable branching ratio into a chargino- state (~ 54% at M2 = 300 GeV). In 
addition, also the channel X3 ~^ X2Z becomes important in this case (with a branching ratio of about 22% at M2 = 
300 GeV). 

B. Branched chain benchmarks 

In the branched chains the decay amplitude is sensitive to the value of mA, while it does not depend on mj. As a 
consequence of this, as far as scenario A is considered, for branched decays we adopt the following benchmark: 
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BR(q^q tt^j) (sequential, benchmark j^ — seq2) 
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FIG. 8: Branching ratios for the sequential process q — > frxi in benchmark ^-seq2 as functions of M2. Each of the first three 
panels corresponds to a different intermediate neutralino state {xi,i = 2, 3, 4) in the sequential chain, while the last panel refers 
to combined branching ratios. The codes for the various curves are as in Fig. |B] 



^-brc : Ml = 10 GcV fi ^ 110 GcV 

tan [3 = 35 tua = 90 GeV ^ = 0, (17) 

while the value of mj is undetermined (for definitcness, in this case we fix mj =150 GeV). This benchmark is depicted 
quahtatively in Figllll and summarized in Table |TT1 The branching ratios for benchmark ^-brc are displayed in 
FigHH where the dotted lines show the branching ratio for the process Xi ~^ Zxi, the short-dashed lines correspond 
to the process Xi ~^ hxi^ the long-dashed lines indicate the process Xi ~^ Hxi, the dot-dashed lines denote the 
process Xi ~* ^Xi and, finally, the solid line indicates the combined branching ratio for the whole decay chain 
q qxi — > <liZ-,h, H, A)xi qbbxi- Note that in this Figure the branching ratios for the process q — > qxi, already 
displayed in Figs. I5I6I for the same set of parameters, are omitted. The combined branching ratio for the decay 
q — > qxi — > qZxi ^e^e Xi, to which only the Z exchange contributes, can be simply calculated by scaling the line 
for Z by a factor 0.22, the ratio of BR{Z ~* e+e") and BR{Z bb) in the Standard Model. 

As shown in Fig[T2l the major contribution to the branching ratio BK{q bbxi) is due to the production of X2, 
which subsequently decays into xi through the production of a Z boson or a A Higgs boson, provided \m2\ is above 
threshold (which implies that M2 > 150 GeV). These two channels have large branching ratios, due to the fact that 
X2 is mainly a Higgsino, and also xi has a sizable Higgsino component. The branching ratio of the channel through 
Z prevails over the branching ratio of the channel through A roughly by a factor 4, since a factor 7-10 due to the 
Lorentz structure of the matrix elements is partially compensated by a factor 0.7-0.5 due to the different coupling 
constants. 
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FIG. 9: 

Figm 



Branching ratios for the sequential process q eexi in benchmark S-seq as functions of M2. Notations are as in 



As far as branched decays are concerned, the most relevant feature of the scenario B is that niA is heavy {niA > 
200 GeV), while tan/3 < 10 is moderate. As a consequence of this, lower branching ratios are expected compared to 
the previous case of benchmark A-hvc. In scenario B we adopt the two following benchmarks: 



5- brcl : tua = 200 GeV 

6- brc2 : tua = 1000 GeV, 



(18) 
(19) 



where, in both cases, fi = —500 GeV, Mi=25 GeV, tan/3=10, jtz^- = 120 GeV and A = 0. These benchmarks are 
depicted qualitatively in FiglTT]and summarized in Table Ull The corresponding branching ratios are shown in FigsflSl 
and [Ml Note that in FiglT4]only the decays Xi ~^ ^Xi ^nd Xi ~^ ^xi are kinematically allowed. The main features 
of these figures may readily be accounted for by arguments similar to the ones previously described. 

C. Expected number of events 

In order to assess the experimental detectability of our scenarios, we now display in Figsll5l- I17l the total production 
cross sections for the process pp — > X ffxi with ^/s = 14 TcV for each of the benchmarks previously introduced (/ = t 
and / = e, for the cases ^-seql, y^-seq2 and S-seq, are displayed in Fig. [TSl and in Fig. [161 respectively; f = b, for 
the cases .4-brc, ;B-brcl and S-brc2, is displayed in Fig. \T7\i . Notice that the case / = e for branched chains can 
simply be derived by scaling the lines of the top-left panel of Fig. [iTl (production of a Z-boson) by a factor 0.22. the 



ratio of BR{Z 



"e") and BR{Z -> bb) in the Standard Model. 
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BR(q^q rr^j) (sequential, benchmark ^ — seq) 
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FIG. 10: Branching ratios for the sequential process q ttxi in benchmark S-seq as functions of M2. Notations are as in 
FigH 



FIG. 11: Quahtative schemes of benchmarks for branched decay chains: ^-brc, i3-brcl, i3-brc2. For each benchmark, extremes 
values for M2 are considered; M2 ~ \n\ and Ah > for scenario A; M2 < M2 ~ and M2 > \fj,\ for scenario B. . 



A rough estimate of the expected events after two year's high-luminosity running at LHC is indicated on the vertical 
axis on the right of these fi gur es. For the calculation of the production cross sections a{pp —^ qq,qq* ,gg,qg) we have 
used the code PROSPINO ^ 

We assume here arbitrarily that the cuts required to extract the SUSY signal from the background will have an 
efficiency of the order of 10%, and a reasonably precise measurement can be performed for a chain for which ^ 100 
events are left after the experimental cuts. 

From the results summarized in Figs llSI - [T7| it turns out that there are very good perspectives for a fruitful 
investigation of the supersymmetric parameter space relevant for light neutralinos at LHC. In particular we notice 
that: 1) The scenario A should be easily explorable both through sequential decay chains and through branched chains. 
2) For the sequential case good perspectives are offered by the ee signal and by the tt one in both benchmarks A- 
seql and ^-seq2. In particular, notice that in the benchmark ^-seql the ee and the fr signals are about of the 
same size, because they take origin from the decay of Xi-i which for M2 > \fJ.\ is dominantly a gaugino; on the 
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BR(q^q bb^^) (branched, benchmark ^ — brc) 
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FIG. 12: Branching ratios for the branched chain in benchmark ^-brc as functions of Af2- Each of the first three panels 
corresponds to a different intermediate neutralino state (Xi)* = 2,3,4) in the branched chain, while the last panel provides 
the branching ratios summed over all Xi's- Notations are as follows. Dotted lines: branching ratio for the process Xt ~^ ^Xii 
short-dashed lines: branching ratio for the process Xi hxi', long-dashed lines: branching ratio for the process Xi ~* Hxi', 
dot-dashed lines: branching ratio for the process Xi ~* ^Xii solid line: combined branching ratios for the whole decay chain 
9 ^ QXi ^ 9(^1 ^1 A)xi — > qbbxi- The branching ratios for the process q —> qxi, already displayed in Figs. I7I8I are omitted 
here. 

contrary, in the benchmark y^-seq2 the f r signal is larger than the ee one, since in this case the process goes mainly 
through intermediate Higgsino-like neutralino states. 3) For the branched case the scenario A gives good measurement 
perspectives for production of a bb pair through Z and h, A Higgs bosons. 4) In scenario B large signals are expected 
in terms of ee and fr pairs in sequential processes; these two signals are comparable in size, since they are generated 
by intermediate neutralino states which are dominantly gauginos. 5) Higher statistics will be required to explore the 
scenario B by branched decay; the most favorable processes are represented by those mediated by Z or h, A Higgs 
bosons when M2 > 

The rates are of course function of the assumed squark mass of 1 TeV. Since the BR for the considered chains do 
not depend on the squark mass, results for different squark masses can be obtained by scaling down the curves by the 
relative squark production cross-section. As the squark mass gets nearer to the gluino mass a significant contribution 
will also come from squarks produced in gluino decays. For instance for a gluino mass of 2 TeV and a squark mass of 
~ 1500(1900) GeV, the scaling factor would be respectively ^ 6(25), and the fraction of events with at least a gluino 
in the initial state would be respectively ^ 18%(55%). 

We stress that here our considerations are simply based on the evaluated total number of events. The actual 
potentiality of investigation at LHC of the present signals will require a detailed analysis in terms of signal to 
background ratios and specific kinematical distributions. This further investigation is beyond the scope of this paper 
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FIG. 13: Branching ratios for the branched chain in benchmark B-brcl as functions of Af2- Notations are as in Fie: ll2l The 
branching ratios for the process q — > qxi, already displayed in Figs. I9I10I are omitted here. 



and will be presented elsewhere [23 1 



CONCLUSIONS 



In this paper we have analyzed the discovery potential of LHC with respect to light neutralinos, i.e. neutralinos 
with a mass < 50 GeV, which arise in supersymmetric models where gaugino masses are not unified at a Grand 
Unified (GUT) scale. These neutralinos have been thoroughly investigated in Refs. [H, under the hypotheses that 
R-parity is conserved and that the lowest neutralino state is the Lightest Supersymmetric Particle (LSP). This LSP 
light neutralino has been proved to be quite interesting as a cold dark matter particle intervening in a number of 
direct and indirect astrophysical effects. 

In particular, in Rcfs. [l|, Q it was derived that the present constraints due to accelerator and other precision 
measurements, together with limits imposed by cosmological observations, concur in confining these neutralinos to 
configurations located in a well delimited part of the supersymmetric parameter space. In other words, the relevant 
region of the SUSY parameter space can be described by a limited number of free parameters, since some of the model 
parameters are essentially frozen by the external constraints. 

Such a situation implies that a few scenarios and relevant conspicuous benchmarks can naturally be singled out. 
This is the strategy that we have used in the present paper in order to explore the discovery potential of LHC as far 
as light neutralinos are concerned. 

The simplicity of the underlying supersymmetric model also allows the derivation of analytic formulae which help 
a lot in understanding the main properties of the spectroscopy of the four neutralino states. The relevant expressions 
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FIG. 14: Branching ratios for the branched chain in benchmark B-brc2 as functions of Af2- Notations are as in Fie: ll2l The 
branching ratios for the process q — > qxi, already displayed in Figs. I9I10I are omitted here. 

have been derived in the first part of the paper. 

Two main scenarios have been introduced: the scenario A, where the stable neutrahno has a mass ^ 10 GeV, 
and the scenario B where 25 GeV (the specification for the values of the other supersymmetric parameters is 

given in Table Within these two scenarios a number of convenient benchmarks have been introduced (see Table 

in. 

Ill the framework of the selected benchmarks we have considered the following decay chains, generated by squarks 
produced in the proton-proton scattering: q qxi qff iffxi (sequential chain), and q —f qxi — > 
q{Z, h, H, A)xi qffxi (branched chain). We limited our discussion to the case in which the gluino is heavier 
than the squark; for definiteness, we have set the SU(3) gaugino mass at the representative value M3 = 2 TeV and 
the squark soft-mass at the value = 1 TeV. Notice that these two parameters are irrelevant in the specification of 
our scenarios inspired by cosmology. 

Branching ratios and number of events expected at LHC have been evaluated for the signals which proved to be 
the most important ones for experimental investigation. 

We arrived at the following main conclusions: 

i) The scenario A should be easily explorable both through sequential decay chains and through branched chains. 
For the sequential case good perspectives arc offered by the ee (or /I/i) signal and by the fr one. For the branched 
case the scenario A gives good measurement perspectives for production of a bb pair through Z and h, A Higgs bosons, 
combined with light lepton pairs through Z. 

ii) In scenario B large signals are expected in terms of ee (or /I/i) and f r pairs in sequential processes. High statistics 
will be required to explore the scenario B by branched decay; the most favorable processes being represented by those 
mediated by Z or h, A Higgs bosons when M2 > |/^|- 
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FIG. 15: Total production cross sections for the processes pp — > Xrfxi in the benchmarks for sequential chains for ^/s = 14 
TeV, as functions of A/2. The decay branching ratios included in the calculation are those displayed in Figs. I6I8I10I The 
notation is as follows. Thin-solid line: ^-seql with mj- = 150 GeV and mediation of fi; thick-solid line: yt-seql with mj- = 150 
GeV and mediation of T2; thin-dotted line: ^-seq2 with mediation of fi; thick-dotted line: ^-seq2 with mediation of f2; 
thin-dashed line: S-seq with mediation of fi; thick-dashed line: S-seq with mediation of T2. On the right vertical axis the 
corresponding total number of expected events is shown, assuming a luminosity of 200 fb^^ (two years of high-luminosity run 
at LHC). 



These results show that LHC has a strong potential in the investigation of the supersymmetric parameter region 
compatible with a light neutralino. Due to the characteristic features of this region, the measurements of LHC should 
easily prove or disprove our model. 

Finally, we wish to recall that our conclusions are based on the evaluated total number of events, only. To ascertain 
the actual potentiality of LHC in the study of our model, the investigation has to be pursued to include analysis of 
the signal/background ratios and of specific kinematical distributions. This further investigation will be presented in 
a subsequent publication [2^. 



19 




FIG. 16: Total production cross sections for the processes pp — > Xeexi in the benchmarks for sequential chains for ^/s — 14 
TeV, as functions of M2 . The decay branching ratios included in the calculation are those displayed in Figs. I5I7I9I The notation 
is as follows. Thin-solid line: ^-seql with mj — 150 GeV and mediation of e^; thick-solid line: ^-seql with m,- = 150 GeV and 
mediation of br; thin-dotted line: ^-seq2 with mediation of bl; thick-dotted line: ^-seq2 with mediation of cr; thin-dashed 
line: S-seq with mediation of cl; thick-dashed line: S-seq with mediation of Sr. On the right vertical axis the corresponding 
total number of expected events is shown, assuming a luminosity of 200 fb~^ (two years of high-luminosity run at LHC). 
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